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^attenuation  were  measured  by  the  radar  system,  with  liquid  water  content  calcu- 
lated from  the  attenuation  information.  The  surface  raingages  were  utilized  to 
obtain  rainfall  rates  and  amounts  while  the  raindrop  spectrometers  yielded 

dropsize  spectra  as  well  as  rainfall  rate  data. 

The  radar  data  indicated  large,  small-scale  variations  for  the 
variables  measured  along  the  radar  radials  of  interest  for  four  rain  periods 
examined.  Both  convective  and  stratiform  precipitation  were  sampled,  the 
convective  rains  having  greater  maxima  of  the  variables  as  well  as  greater  time 
and  space  variability  for  the  measured  parameters. 
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LINEAR  PRECIPITATION  CHARACTERISTICS  IN  THE  ATMOSPHERE 

INTRODUCTION 

This  AFGL  sponsored  contract  was  designed  to  investigate  the  relationship 
between  precipitation  cha rac ter i s i t i cs  along  ray  paths  in  the  atmosphere  and 
those  along  the  ground. 

The  observational  program  associated  with  the  research  utilized  a dual 
wavelength  (3  and  10  cm)  radar  system,  raindrop  spectrometers  and  recording 
raingages  to  provide  data  for  defining  variations  in  liquid  water  content,  radar 
reflectivity  and  surface  rainfall  characteristics.  Potential  field  operational 
periods  during  1 9 7 3 - 7 5 were  in  spring  and  fall  in  eastcentral  Illinois  and  in 
the  St.  Louis,  Missouri  region  in  summer.  The  primary  goal  was  to  collect  a 
sufficiently  large  quantity  of  data  from  varying  precipitation  regimes  to  permit 
the  definition  of  statistical  distributions  of  the  above  parameters  along  radar 


radials  in  nearly  instantaneous  time  intervals,  and  the  development  of  relationships 
between  these  distributions.  Data  collection  was  often  limited  by  the  use  of  much 
of  the  needed  equipment  (radar  and  raingages)  for  other  grants  and  contracts 
which  had  higher  operational  priorities. 


OBSERVATIONAL  AND  ANALYTICAL  PROGRAM 

I nt  roduc  t i on 

The  primary  purpose  for  development  of  the  dual  wavelength  radar  was 
as  a hail  detection  device,  and  its  development  was  funded  by  a separate, 
joint  NSF  grant  to  the  Water  Survey  and  the  University  of  Chicago.  This  radar 
is  referred  to  as  the  CHILL  radar  (CH-Chicago,  I LL- I 1 1 i no i s) . Its  presence  in 
Illinois  therefore  was  to  support  hail  research  with  ground  truth  from  a surface 
network,  of  recording  raingages  and  hail  pads.  As  a result,  the  use  of  the  radar 
for  this  AFCRL  effort  was  contingent  on  the  radars  availability  as  established 
by  the  NSF  grant  which  financially  supported  radar  operations. 


Field  operations  in  197^  and  1975  related  to  the  line-oriented 
recording  raingages  and  raindrop  spectrometers  were  funded  by  this  contract 
and  were  operated  in  conjunction  with  the  radar  system. 

Data  Acquisition  Systems 

Radar . Since  variations  in  the  liquid  water  content  (LWC)  are  of 
prime  interest  to  this  AFGL  study,  the  dual  wavelength  radar  utilized  (when 
available)  provided  an  adequate  instrumentation  system  since  it  can  measure 
attenuation  rate  (as  well  as  radar  reflectivity);  this  parameter  is  necessary 
for  the  calculation  of  LWC.  This  radar  provides  information  concerning  the 
gross  features  of  microwave  attenuation  and  attenuation  statistics. 

The  detection  of  hail  and  the  determination  of  hail  size  was  the  initial 
impetus  for  the  development  of  this  dual  wavelength  radar  system.  Theoretically, 
hail  detection  is  accomplished  with  the  dual  wavelength  radar  by  measuring  the 
average  echo  power  returned  at  10  cm  and  3 cm  wavelengths  from  matched  beams 
illuminating  volumes  containing  cloud  and  precipitation  particles.  The  range 
derivative  of  the  logarithm  of  the  ratio  of  the  average  echo  powers  at  the 
10  cm  and  3 cm  wavelengths  is  then  calculated  and  used  as  the  hail  identifying 
signal.  Determination  of  attenuation  results  from  the  above  procedure  and 
attenuation  rate  may  then  be  used  to  obtain  LWC  with  the  aid  of  an  empirically 
determined  equation  relating  the  two  parameters. 

Previous  methods  of  LWC  determinations  utilizing  radar  involved  the 
measurement  of  returned  power  from  a single  wavelength  radar  from  which  the 
radar  reflectivity  factor  (Z)  was  calculated  using  the  radar  equation; 
attenuation  rate  was  then  calculated  from  empirical  relationships.  Additional  errors 
result  from  utilizing  this  method  since  attenuation  is  empirically  determined. 
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The  dual  wavelength  technique  measures  attenuation.  Further  details  concerning 
LWC  determination  utilizing  dual  wavelength  radar  may  be  found  in  a paper  by 
Eccles  and  Mueller  (1971). 

Since  the  radar  system  was  originally  designed  and  utilized  as  a tool 
for  hail  detection  other  federally  sponsored  studies,  the  field  periods 
for  this  AFGL  study  were  scheduled  to  coincide,  when  possible,  with  the 
availability  of  the  radar  in  east  central  Illinois  when  it  was  not  being 
utilized  elsewhere.  The  funding  for  operation  of  the  radar  was  provided  by 
separate  qrants  (NSF  DES75- 1 422 1 and  NSF  GK-37859) . Reliability  of  the  system 
was  not  as  high  as  desirable.  However,  it  must  be  realized  that  this  is  not 
an  "off  the  shelf"  item  and  still  had  some  inherent  problems  often  associated 
with  newly  developed  instrumentation.  Most  of  the  difficulties  were 
associated  with  signal  processing  and  recording  of  the  radar  information.  The 
data  set  obtained  in  spring  1975  was  more  usable  than  that  from  previous 
periods  and  provided  the  bulk  of  the  radar  information  used  in  analyses 
reported  on  further  in  this  report. 

Digitized  radar  data  were  recorded  on  magentic  tape  when  rain  was 
occurrring  within  the  surface  instrumentation  network.  In  addition  to  tape 
recording  the  radar  information,  a remote  PP I scope (p 1 an-pos i t i on  indicator) 
was  photographed  at  approximately  l-minute  intervals  as  a backup  and  as  a 
check  for  the  tape  recording  system.  The  radar  antenna  scanning  was  programmed 
to  accomodate  the  hail  detection  work  as  a first  priority.  Antenna  elevation 
angle  generally  was  incremented  by  approximately  1°  at  the  end  of  each  360° 
sweep  to  a maximum  elevation  determined  by  the  echo  tops  and  the  distance 
between  the  radar  and  the  echoes  of  interest.  Time  between  sweeps  for  a 


given  elevation  angle  was  variable  as  a result  of  the  above  antenna  programming 
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mode.  However,  this  time  interval  generally  did  not  exceed  5 minutes. 

Data  collection  continued  until  the  precipitation  left  the  network. 

Raindrop  Spectrometers.  Six  raindrop  spectrometers  were  utilized  in 
the  data  collection  phase  to  obtain  raindrop  size  distributions  at  the 
surface  associated  with  the  rains  of  interest.  This  instrument  was  designed 
and  developed  at  the  Water  Survey  (with  ERDA  funding)  over  the  past  five 
years  with  a number  of  modifications  being  made  during  this  time.  The 
raindrop  spectrometer  is  a device  which  senses  the  momentum  of  raindrops, 
assumed  to  be  falling  at  terminal  velocity  at  ground  level,  and  converts 
these  signals  to  digits  which  are  proportional  to  drop  size.  A transducer 
in  the  system  transforms  the  mechanical  impulse  generated  by  a drop  striking 
the  sensor  head  to  an  electronic  signal  which  is  then  digitized  and  recorded 
on  magnetic  cassette  tape  for  further  processing.  Instrument  reliability 
has  slowly  improved,  however,  as  was  indicated  in  previous  quarterly  reports, 
operational  difficulties  were  encountered  during  the  data  collection  periods 
of  this  contract.  As  a result,  the  amount  of  usable  information  is  small. 

Since  two  of  the  parameters  calculated  from  drop  spectra  are  radar 
reflectivity  factor  (Z)  and  liquid  water  content,  these  values  measured  at 
the  ground  may  be  compared  to  the  above  ground,  radar  measured  Z and  the  LWC 
calculated  from  the  radar  measured  attenuation.  Rainfall  rate  is  also 
calculated  from  the  drop  spectra,  usually  over  60  second  intervals. 

Recording  Raingages.  Six  weighing  bucket  type  recording  raingages  were 
used  in  addition  to  the  raindrop  spectrometers  to  yield  rainfall  rates  and 
amounts  at  the  surface.  Operation  of  the  gages  in  Illinois  was  funded  by 
this  contract.  The  clock  mechanism  in  each  was  geared  so  that  rates  could 
be  measured  to  I-minute  time  intervals.  The  gage  also  was  fitted  with 
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a 32  cm  diameter  top  as  opposed  to  the  standard  20  cm  top  to  improve  amount 
resolution.  The  accuracy  of  these  short-time  rate  measurements  is  limited 
by  the  mechanical  and  hydraulic  response  limitations  of  the  instrument  as  well 
as  by  the  human  error  involved  in  digitizing  the  raingage  chart  information. 

The  gage  information  is  recorded  on  a 29  cm  long  chart  with  a clock  cylinder 
geared  to  rotate  once  each  6 hours.  On  these  charts,  1 minute  is  equivalent 
to  0.8  mm  along  the  horizontal  dimension  and  2.5  mm  vertically  is  equivalent 
to  1.0  mm  of  precipitation.  These  resolution  problems  do  not  exist  with  the 
rates  determined  from  raindrop  spectra. 

Field  Operat i ons 

The  spring  and  fall  field  operational  periods  (Spring  197^  and  1975. 

Fall  1973  and  197^)  typically  consisted  of  8-10  weeks  of  activity  involving 
personnel  associated  with  the  radar-hail  detection  program  as  well  as  those 
connected  with  this  effort.  Two  individuals  were  scheduled  for  duty  on  each 
day  during  the  period  and  were  responsible  for  examining  the  weather  information 
as  provided  by  a weather  facsimile  machine  and  c i v i 1 i an-a i rways  circuit  teletype 
printer.  The  synoptic  conditions  were  monitored  throughout  the  day  so  that 
potential  precipitation  situations  would  be  recognized  and  necessary  preparations 
made  regarding  support  personnel  and  equipment.  When  precipitation  over  the 
ra i ngage-ha i 1 pad-spect rome ter  network  was  imminent,  the  radar  was  activated. 

The  spectrometers  and  gages  were  operational  at  all  times  and  were  routinely 
visited  by  field  technicians  at  intervals  which  assured  sufficient  servicing 
for  the  desired  time  resolution,  and  also  insured  that  "missed"  precipitation 
occurrences  were  kept  to  a minimum.  These  sites  were  serviced  as  soon  after 


a rain  as  possible. 


-8- 


Data  Collection  Procedure  and  Format 

The  instrumentation  network  for  this  project  in  east  central  Illinois 
was  activated  in  conjunction  with  one  designed  to  test  the  dual  wavelength 
radar  system  as  a tool  for  hail  detection.  The  network  shown  in  Figure  1 
depicts  the  locations  of  the  six  raindrop  spectrometers  and  located  raingages 
with  respect  to  the  radar  and  to  the  surface  hail  network.  The  closest 
spectrometer-raingage  site  to  the  radar  was  approximately  35  Km  with  the 
farthest  one  being  approximately  58  Km. 

There  were  four  data  collection  periods  related  to  this  3~year  effort 
with  two  occurring  in  spring  and  two  during  fall.  The  most  recent  one  in 
Spring  1975  was  utilized  exclusively  for  the  data  analysis  phase  since  data 
quality  and  quantity  were  significantly  better  than  in  previous  periods. 

The  data  from  four  rain  events  during  the  Spring  75  period  will  be  presented. 
Operational  problems  with  the  radar  and  spectrometers  greatly  limited  the 
value  of  the  1973  and  197^,  and  portions  of  the  1975  data. 

Radar  data  were  recorded  on  magnetic  tape  with  Z and  A (attenuation) 
values  measured  for  each  150  m-)ength  range  bin  along  radials  of  interest 
emanating  from  the  radar  location.  Data  were  recorded  for  radials  (azimuths) 
270°,  271°  and  272°  (see  Fig.  1),  the  raingages  and  sepct romete rs  being  along 
the  2 7 1 ° radial.  These  recorded  data  included  0.5°,  1 • 5° , 2.5°,  3-5°  and  k° 
elevation  angles  along  each  radial;  height  of  the  radar  beam  above  ground  of 
these  elevation  angles  as  a function  of  range  is  shown  in  Figure  2.  These 
elevations  allow  the  determination  of  LWC  for  heights  up  to  the  freezing 
level.  Averaging  of  Z and  A values  was  performed  over  two  and  three  consecutive 
range  bin  intervals  along  the  radials  in  order  to  obtain  representative 
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samples.  The  LWC  values  were  obtained  from  attenuation  rate  through  the 
empirical  equation: 

M = 2.23A^’^^  (Eccles  and  Mueller,  1971) 

where  M is  LWC  in  g/m^  and  A is  attenuation  rate  in  db/Km.  This  equation 
was  derived  from  LWC  and  A values  calculated  from  approximately  A00  independent 
cubic  meter  samples  of  raindrop  data  obtained  with  a photographic  technique. 

There  are  limitations  to  the  use  of  this  expression  which  are  important 
to  note.  Erroneously  high  LWC  values  may  result  when  this  expression  is  used 
for  radar  volumes  which  contain  solid  (ice)  precipitation  particles.  Therefore 
caution  should  be  used  when  interpreting  LWC  data  from  convective  systems 
that  may  contain  hail  or  when  t fie  radar  beam  is  above  the  freezing  level. 
However,  in  the  absence  of  solid  precipitation,  a good  estimate  of  LWC  may 
be  made,  particularly  below  cloud  base  in  a rain-only  environment  as  well  as 
in-cloud  in  the  presence  of  cloud  droplets  and  rain. 

It  is  of  value  to  know  what  can  be  expected  as  an  upper  limit  to  LWC 
in  clouds  and  precipitation  in  order  to  have  an  empirical  reference  with 
which  to  compare  data  obtained  using  the  above  scheme.  Roys  and  Kessler  (1966) 
indicated  a value  of  approximately  A3  g/m  as  a maximum  in-cloud  measurement 
of  LWC  during  cloud  penetrations  into  thunderstorms  with  an  F-100F  aircraft. 
Surface  values  of  LWC  calculated  from  raindrop  spectra  obtained  with  a 
photographic  technique  approach  a maximum  of  approximately  29  g/m^  tor 
spectra  recorded  with  the  raindrop  camera. 

Surface  rainfall  rates  and  amounts  were  obtained  from  weighing  bucket 
type  recording  raingages;  these  gages  are  the  type  used  widely  by  the  National 


Weather  Service.  Rate-time  resolution  may  be  varied  by  changing  the  qear 
ratio  controlling  the  rotation  speed  of  the  chart  on  which  the  information  is 
recorded.  This  permits  a maximum  resolution  of  1-minute  for  rate  determinations 
when  the  proper  gear  ratio  is  utilized.  However,  as  indicated  previously,  the 
accuracy  of  these  short-time  rates  with  this  rainqage  system  is  affected  by 
desiqn  and  human  error  factors  which  permit  significantly  more  error  than 
1-minute  rates  determined  from  raindrop  spectra.  As  may  be  seen  in  Figure  I, 
six  recording  raingaqes,  each  spaced  approximately  5 km  apart  along  the  271° 
radial  from  the  radar,  were  utilized  for  this  study. 

The  raingage  charts  pertaining  to  the  periods  of  interest  were  digitized 
with  an  Autotrol  chart  reader.  A computer  program  was  later  used  to  convert 
these  digitized  values  to  1-minute  rainfall  rates  and  amounts.  This  information 
is  then  related  to  the  corre -.pond  i ng  radar  data. 

Raindrop  data  were  recorded  on  cassette  tape  for  subsequent  data  reduction 
by  computer.  In  addition  to  recording  raindrop  siqnals  on  tape,  date  and  time 
were  also  recorded  from  a built-in  electronic  clock;  this  allows  "rain  time" 
to  be  determined.  The  system  tape  recorder  is  activated  during  rain  only 
althouqh  the  clock  runs  continuously.  Considerable  difficulty  was  encountered 
with  spectrometer  operations  during  the  data  collection  periods.  Operations 
in  the  Spring  1975  period  yielded  the  most  usable  drop  spectra  data. 

Data  Reduction  and  Ana  lysis 

There  were  four  rain  periods  during  Spring  1975  from  which  adequate 
radar  and  rainqage  information  were  obtained.  Spectrometer  data  collected 
during  this  time  were  sporadic  due  to  a number  of  instrument  malfunctions. 

A description  of  each  period  follows. 
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18  April.  Approximately  35  minutes  of  radar  data  were  collected  as 
light  precipitation  associated  with  an  approaching  cold  front  intersected 
the  271°  radial.  Rain  amounts  recorded  by  the  surface  raingages  during 
this  time  were  light  with  measurable  precipitation  at  3 of  the  6 gages  (0.8( 

0.5  and  0.5  mm).  There  were  no  indications  that  this  rain  was  associated 
with  thundershowers  although  a previous  rain  occurring  approximately  2 hours 
earlier  did  contain  thundershowers.  The  freezing  level  just  prior  to  the 
rain  was  at  3.3  Km.  From  Figure  2,  elevation  angles  >2°  would  contain  radar 
volumes  above  the  freezing  level  in  the  western  portion  of  the  network. 

23  April.  Light  precipitation  associated  with  an  approaching,  slow 
moving  cold  front  occurred  over  the  raingage-spectrometer  network  during  the 
forenoon.  The  rain  amounts  were  fairly  uniform  across  the  line  with  all 
gages  receiving  rain  amounts  between  3.0  and  4.6  mm.  There  were  no  indications 
than  any  of  the  precipitation  over  the  line  was  convective  in  nature.  The 
freezing  level  prior  to  the  rain  was  at  approximately  3.7  Km.  One-minute 
rainfall  rates  did  not  exceed  14  mm/hr  along  the  line  with  average  rates 
between  5-8  mm/hr. 

30  April.  Radar  data  were  recorded  during  a 65  minute  period  on  this 
day.  The  rainfall,  associated  with  an  active  cold  front,  was  convective  in 
nature  with  thundershowers,  and  ail  gages  along  the  radial  recorded 
precipitation.  Rain  amounts  varied  from  11.9  mm  at  one  end  of  the  raidal 
decreasing  to  2.0  mm  at  the  other.  Several  1 -minute  rates  exceeded  50  mm/hr 
with  a maximum  rate  during  the  period  of  159  mm/hr.  The  freezing  level  during 
precipitation  was  approximately  2.9  Km.  For  a more  detailed  analysis  of  this 


case  see  Appendix  D. 
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30  Hay.  Approximately  85  minutes  of  radar  data  were  recorded  during 
this  strongly  convective  situation  associated  with  a rapidly  moving  cold 
front.  Rain  amounts  along  the  271°  radial  ranged  from  7-6  mm  to  11.4  mm 
during  the  above  period  with  1-minute  rates  as  high  as  60  mm/hr  observed.  Hail 
also  was  recorded  at  the  ground  at  3 of  the  6 rai ngage-spectrometer  sites 
during  the  rain  although  the  number  density  of  stones  was  small.  The  freezing 
level  prior  to  the  rain  was  at  approximately  3-7  km. 

Average  values  for  LWC  and  Z are  presented  in  Appendix  A for  18  April, 

23  April,  30  April  and  30  May  1975-  The  data  shown  are  for  all  radar  radials 
from  270°  to  212 0 for  elevation  angles  0.5°  to  4.5°;  2 and  3 - b i n moving  averages 
are  both  shown.  Average  values  for  R as  determined  from  the  six  raingages 
located  along  the  271°  azimuth,  are  shown  in  Appendix  B for  times  concurrent 
with  the  radar  data.  Average  raindrop  spectra  when  available  are  presented  in 
Appendix  C for  the  above  periods. 

The  Z,  A,  and  LWC  values  along  the  radials  of  interest  were  averaged 
over  2 and  3"bin  intervals  as  a moving  average  along  each  radial.  Each  value 
then  represented  a radar  volume  300  or  450  m in  length  by  1°  in  beam  width. 

The  number  of  bins  to  be  averaged  for  a representative  sample  was  not  entirely 
clear  initially.  When  2 and  3"bin  radar-determined  LWC  values  are  compared 
for  the  two  convective  days,  (Fig.  3),  average  maximum  values  are  greater  for  the 
2-bin  data.  The  2-bin  data  also  show  more  variance  in  the  average  maxima.  However, 
examination  of  Z shows  no  observable  difference  in  the  2 and  3_bin  results. 

The  significant  difference  in  the  2 and  3"bin  LWC  averages  results  using 
a slightly  different  data  set  for  each  determination.  A greater  number  of  bin 
averages  are  discarded  for  the  3'bin  averaging  than  2-bin  averaging  because  bins 
containing  undefined  reflectivity  values  are  not  used  in  the  LWC  calculations. 

Since  2-bin  averaging  yields  better  space  resolution,  this  approach  is  considered 


more  des i reab I e . 


Another  aspect  of  LWC  variations  may  be  examined  with  the  l-minute  rainfall 
rate  data  obtained  with  the  recording  raingages.  An  excellent  correlation 
exists  between  rainfall  rate  and  liquid  water  content  calculated  from  surface 
determined  raindrop  spectra  of  one  minute  durations.  Mueller  and  Sims  ( 1964) 
found  that  LWC  and  l-minute  rainfall  rate  are  highly  correlated  when  these  values 
are  calculated  from  raindrop  distributions  obtained  from  a photographic  device 
(Mueller  and  Sims,  1966).  The  rel at i onsh i p between  these  variables  was  determined 
for  several  locations  around  the  world,  and  the  relationship  changed  insignificantly 
from  location  to  location.  When  this  relationship  is  determined  from  l-minute 
spectra  obtained  with  the  raindrop  spectrometer,  it  too  is  not  significantly 
different.  The  expression  relating  LWC  and  R is: 

LWC  = .042  R (1) 

with  a standard  error  of  estimate  (SE)  on  LWC  generally  <0.1;  a plot  of  LWC  versus 
R for  data  obtained  with  the  raindrop  camera  and  drop  spectrometer  is  shown 
in  Figure  4. 

When  the  expressions  used  to  calculate  LWC  and  R are  examined,  it  may  be 
seen  that  the  terminal  velocity  of  the  raindrops  in  a drop  distribution  is  the 
only  variable  that  is  not  common  to  both  expressions. 

Rainfall  rate  as  calculated  from  spectrometer  data  is  given  by 

D=6 . 0 

R(mm/hr)  = K £ IT  D^Nq 
D=0.1  Z 

where  K is  a constant  which  depends  on  the  sensor  area  and  sample  time,  and 

D=6 . 0 HD^Nq 
LWC  (g/rrr)  = T.  6 


where  v-j-pis  the  terminal  velocity  of  drop  size  D. 
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Therefore,  the  variability  between  spectra  and  the  associated  terminal 
velocities  is  responsible  for  the  variance  in  the  LWC , R relationships. 

However,  since  the  median  terminal  velocity  from  spectra  to  spectra  varies 
little,  the  SE  is  small.  As  a result,  LWC  values  between  cloud  base  and 
ground  can  be  estimated  with  reasonable  accuracy  by  use  of  surface  determined 
rainfall  rate  in  equation  (1). 

This  approach  can  also  serve  as  a check  for  LWC  values  obtained  with  the 
dual  wavelength  technique  since  rainfall  rate  was  recorded  at  the  six  locations 
indicated  in  Figure  1.  It  must  be  realized  that  a comparison  of  this  type 
necessitates  averaging  a number  of  values  in  space  and  time,  since  radar 
parameters  pertaining  to  a volume  are  compared  to  surface,  point  measurements 
of  rainfall  rate.  A plot  of  LWC  as  calculated  form  1-minute  raindrop  spectra, 
radar-determined  LWC  and  LWC  calculated  from  equation  (1)  using  the  l-minute 
raingage  determined  rainfall  rates,  is  shown  in  Figure  5a.  These  data  relate 
to  information  pertaining  to  raingage  location  67  (Fig.  1)  for  the  30  April  1975 
rain  period.  Values  of  radar-determined  LWC  are  generally  higher  than 
corresponding  surface  values  for  this  data  set.  This  also  occurred  at  location 
65  shown  in  Figure  5b  although  there  is  a wide  discrepancy  during  a portion 
of  the  rain  between  rainfall  rate  (and  resulting  LWC)  from  the  raingage  versus 
the  spectrometer.  An  electronic  saturation  problem  in  this  spectrometer  did  not 
allow  the  recording  of  rain  rates  exceeding  approximately  60  mm/hr.  Average 
raindrop  spectra  for  the  above  rain  periods  are  shown  in  Figure  6.  The  LWC 
discrepancy  between  radar  and  surface-determined  LWC  likely  resulted  from  the 
large  difference  in  sampling  volume  as  well  as  the  difficulty  in  associating,  in 
time  and  space,  a radar  volume  aloft  with  point  measurements  on  the  ground.  The 

obvious  advantage  to  radar-measured  LWC  as  is  also  the  case  with  rainfall  estimated 
from  radar  is  the  large  areal  coverage  with  radar  as  opposed  to  point  measurements 


obtained  from  surface  instrumentation. 
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SUMMARY 

Liquid  water  content  measurements  obtained  from  dual  wavelength  radar 
data  indicate  that  the  averaging  technique  used  is  important  to  the  manner 
in  which  the  data  are  interpreted.  In  a study  such  as  this,  aimed  at  learning 

what  are  the  small  scale  variations  in  LWC  and  Z along  lines,  the  limiting 
factor  becomes  the  smallest  volume  size  which  can  be  examined  for  the  vari- 
ables of  interest.  Both  2-bin  and  3-bin  averaging  of  the  radar-determined 
average  maximum  LWC  are  presented  in  the  appendices  for  the  four  rainfall 
events  sampled.  Regarding  the  more  appropriate  averaging  technique  for  Z,  there 
is  no  significant  difference  in  the  average  2-bin  and  3_bin  maxima  or  standard 
deviation  for  the  rain  periods  analyzed.  This  indicates  that  the  smaller 
total  volume  (2-bin  averaging)  is  probably  large  enough  for  obtaining  represent- 
ative values  of  Z. 

A technique  for  estimating  LWC  from  surface  point  rainfall  rates  was 
also  presented  based  on  raindrop  distributions.  Since  wide  variations  exist 
in  short  interval  rainfall  rates  and  drop  spectra,  the  resulting  LWC  values 
also  show  this  tendency. 

Examination  of  the  data  sets  reveals  that,  as  expected,  there  are  wide 
fluctuations  in  LWC  and  Z as  a function  of  range  along  lines  (radar  radials). 
These  variations  are  not  as  extensive  for  more  stratiform  rains  as  the  18  April 
and  23  April  cases  which  had  no  apparent  convective  elements.  Average  maxima 
as  well  as  variance  increase  for  LWC  and  Z as  the  precipitation  becomes  more 
convective  in  nature. 

Although  hail  was  observed  on  the  ground  during  the  30  May  rain,  the 
LWC  values  are  in  line  with  the  30  April  convective  rain  when  no  hail  was 
recorded.  Apparently,  the  number  of  hailstones  per  unit  volume  was  sufficiently 


small  to  be  "averaged  out"  in  the  statistics.  It  is  also  possible  that  the 
hailfall  was  of  a very  short  duration  and  was  not  occurring  during  the  recorded 
radar  sweeps.  Since  there  was  no  time  information  concerning  the  hail  data, 
this  cannot  be  verified. 

The  major  goal  of  this  study  was  to  obtain  a large  data  sample  of  radar 
and  surface  rainfall  information.  It  is  apparent  from  the  limited  amount  of 
data  reported  on  herein  that  this  goal  was  not  fully  realized.  A major  limit- 
ation arose  from  the  limited  amount  of  funding  allotted  for  the  effort,  in- 
adequate to  support  extensive  radar  operations.  Radar  data  collection  was 
possible  only  when  the  dual  wavelength  radar  system  was  operating  in  east 
central  Illinois  as  dictated  by  funding  for  a separate  project.  There  was 
no  AFGL  funding  allocated  specifically  for  radar  operations;  therefore,  opera- 
tional priorities  were  weighted  towards  the  NSF-funded  radar  effort.  The 
AFGL  radar  effort  was  in  essence  "piggy-backed"  onto  the  other  project.  In 
addition,  poor  quality  of  much  of  the  radar  data  obtained  was  a contributing 
factor  to  the  very  limited  data  bank.  Difficulties  also  with  raindrop  spectro- 
meter operations  contributed  to  the  poor  quality  and  quantity  of  surface 
drop  spectra  information. 

Obviously,  firm  conclusions  cannot  be  drawn  from  the  small  amount  of 
data  collected  and  analyzed.  In  the  event  information  of  this  nature  is 
required  in  the  future,  a similar  effort  would  be  more  successful  if  funds  for 
radar  operations  were  provided. 

Many  of  the  operational  difficulties  encountered  with  the  dual  wave- 
length radar  and  the  raindrop  spectrometers  resulted  from  problems  often 
associated  with  newly  developed  instrumentation.  Most  of  these  have  been 


corrected . 
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Anpendix  A 
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Ti  me 

Mean 

Rainfall  Rate 
(mm/h  r) 

Mini  mum 
Rate 

Max i mum 
Rate 

No . of 
Points 

1 1 15 

2.66 

0.51 

5-73 

005 

1 1 16 

3.65 

0.51 

6.01 

005 

1117 

5-42 

0.56 

13.62 

005 

1118 

5-25 

0.56 

11.59 

005 

1119 

4.13 

1-27 

6.99 

005 

1120 

4.57 

2.99 

6.92 

005 

1121 

4.66 

3-23 

5-65 

005 

1122 

5.37 

3-23 

7-13 

005 

1123 

5-50 

2.56 

8.00 

005 

1124 

5-58 

2.26 

7.42 

005 

1125 

5-38 

2.26 

6.49 

005 

1126 

5.15 

3.78 

6. 10 

005 

1 127 

4.43 

1.73 

6.87 

005 

1128 

4.37 

1.73 

7-53 

005 

1129 

4.69 

1-73 

6.95 

005 

1130 

4.59 

1 -73 

6.29 

005 

1131 

4.27 

1.73 

6.29 

005 

1132 

3.78 

1 .89 

6.29 

005 

1133 

4.43 

2.02 

7.02 

005 

1134 

4.63 

2.02 

8.06 

005 

1135 

4.67 

2.02 

11.34 

005 

1136 

2.19 

1 .29 

3-21 

005 

1137 

2.29 

1 .29 

3.69 

005 

1138 

2.37 

1.85 

3.69 

005 

1139 

2.47 

1 .42 

4.97 

005 

1 140 

4. 12 

0 . 66 

7-99 

005 

1141 

3.65 

0 . 66 

7.47 

005 

1142 

3.17 

0.59 

7.64 

005 

1143 

2.79 

0.24 

7.70 

005 

1 144 

2.90 

0.24 

7-70 

005 

1145 

4.17 

1.42 

7-32 

005 

1146 

3-05 

1.42 

7-25 

005 

1147 

2.50 

0.60 

7-25 

005 

1148 

1.98 

0.00 

5.75 

005 

1149 

1.40 

0.00 

4.69 

005 

1150 

1.66 

0.00 

3-70 

005 

k A 
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30  Apri ] 1975 


Mean 


Rainfall  Rate 

Mini  mum 

Max i mum 

No . of 

T i me 

(mm/h  r) 

Rate 

Rate 

Points 

15^5 

12.57 

12.57 

12.57 

001 

1546 

5-96 

5.96 

5-96 

00  1 

1547 

3-56 

3-56 

3-56 

001 

1548 

3.69 

3-69 

3-69 

001 

1549 

6.19 

6.19 

6 . 19 

001 

1550 

10.57 

0.24 

20.90 

002 

1551 

20.03 

0.00 

40.07 

002 

1552 

23.81 

0.00 

47.63 

002 

1553 

16.39 

0.00 

32.78 

002 

1554 

20.83 

0.00 

41  .66 

002 

1555 

23.80 

0.00 

47-61 

002 

1556 

27.26 

1.48 

53-05 

002 

1557 

44.9) 

2.49 

87.33 

002 

1558 

80.54 

2-34 

158.74 

002 

1559 

49.06 

0.82 

97.30 

002 

1600 

9-94 

3-37 

16.52 

002 

1601 

9.48 

4.84 

14.  12 

002 

1602 

5.96 

5.  13 

6.80 

002 

1603 

10.89 

5.71 

16.07 

002 

1604 

9-94 

0.00 

23.76 

003 

1605 

14.55 

0.22 

37-37 

003 

1606 

19-93 

1.31 

55-53 

003 

1607 

12.37 

0.00 

28.35 

004 

1608 

7-75 

0.00 

26.42 

005 

1609 

10.70 

0.00 

51-99 

005 

1610 

13-36 

0.00 

66.16 

005 

1611 

13.92 

0.00 

69-27 

005 

1612 

16.07 

0.00 

53.63 

005 

1613 

17.69 

0.00 

49-42 

005 

1614 

13-92 

0.00 

40.79 

005 

1615 

11.32 

0.00 

41  .49 

005 

1616 

9-48 

0.00 

36.67 

005 

1617 

5.75 

0.00 

26 . 86 

005 

1618 

7-34 

0.00 

32.85 

005 

1619 

12.57 

0.00 

52.23 

005 

1620 

10.51 

0.00 

41.61 

006 

1621 

10.45 

0.00 

41  .06 

006 

1622 

8.51 

0.00 

21  . 10 

006 

1623 

5.68 

0.00 

26.13 

006 

1624 

5.06 

0.00 

21  .81 

006 

1625 

4.53 

0.00 

20.20 

006 

1626 

7. 06 

0.00 

34.24 

006 

1627 

10.82 

0.  18 

55.60 

006 

1628 

7-32 

0.18 

32.37 

006 

1629 

1.51 

0.18 

3-78 

006 

1630 

I.89 

0. 18 

3-78 

006 

1631 

2.17 

0.18 

3-99 

006 

T i me 

632 

633 
163^ 

1635 

1636 

1637 

1638 

1639 

1640 

1641 

1642 

1643 

1644 

1645 

1646 

1647 

1648 

1649 

1650 


Mean 

Ra i nfa 11  Rate 
(mm/h  r) 

Mini  mum 
Rate 

Maxi  mum 
Rate 

No . of 
Points 

3-07 

0.18 

8.20 

006 

5.03 

0.18 

17.01 

006 

14.03 

0.18 

65-17 

006 

9-35 

0. 18 

39-27 

006 

3-39 

0.18 

6.83 

006 

2.  1 1 

0.  18 

6.59 

006 

1.47 

0.00 

7.21 

006 

1.30 

0.00 

6.43 

006 

0.93 

0.00 

2.67 

006 

0.96 

0.  18 

2.67 

006 

0.96 

0.  18 

2.67 

006 

0.96 

0. 18 

2.75 

006 

1 . 10 

0.18 

2.92 

005 

0.96 

0.00 

2.92 

005 

3.80 

0.18 

13.35 

004 

7- 06 

0.00 

26.29 

004 

4.53 

0.00 

16.  19 

004 

3.48 

0.4S 

10 .51 

004 

1.77 

0.45 

3-30 

003 

T i me 


505 

506 

507 

508 

1509 

1510 

1511 

1512 

1513 

1514 

1515 

1516 

1517 

1518 

1519 

1520 

1521 

1522 

1523 

1524 

1525 

1526 

1527 

1528 

1529 

1530 

1531 

1532 

1533 

1534 

1535 

1536 

1537 

1538 

1539 

1540 

1541 

1542 

1543 

1544 

1545 

1546 

1547 

1548 

1549 

1550 


Mean 

Rai nfal 1 Rate 
(mm/hr) 

Mini  mum 
Rate 

Max i mum 
Rate 

No . of 
Poi nts 

0.87 

0.87 

O.87 

001 

1.80 

1 .80 

1 .80 

001 

37.78 

37-78 

37.78 

001 

46.65 

46.65 

46.65 

001 

30.  18 

30.18 

30.  18 

001 

13.02 

13-02 

13.02 

001 

10.38 

10.38 

10.38 

001 

7.66 

7 . 66 

7-66 

001 

7 . 66 

7.66 

7.66 

00  1 

7.94 

7-94 

7.94 

001 

8.30 

8.30 

8.30 

001 

9-07 

9-07 

9-07 

001 

1 1 .03 

1 1 .03 

1 1 .03 

001 

1.44 

0.00 

2.88 

002 

0.93 

0.00 

1.87 

002 

3.20 

0.50 

5-91 

002 

21  .08 

12.84 

29-33 

002 

22.22 

7-99 

36.45 

002 

5.89 

0.74 

1 1 .05 

002 

2.59 

0.00 

5-19 

002 

0.16 

0.00 

0-33 

002 

0.11 

0.00 

0.22 

002 

4.31 

0.22 

10.99 

003 

1 1.06 

0.22 

31.21 

003 

5.59 

0.00 

15.01 

003 

4.04 

1 .29 

7.21 

003 

3.48 

2.96 

3.84 

003 

4.53 

2.45 

7-31 

003 

6.55 

1.48 

14.86 

003 

6.57 

1.43 

15-18 

003 

4.97 

1-37 

10.44 

003 

5.73 

1.37 

12.65 

003 

9.69 

1.37 

19-59 

00  4 

10.53 

1.37 

34.32 

004 

8.51 

1 . 16 

27.89 

00  4 

4.51 

0.34 

13.75 

004 

2.01 

0.00 

4.04 

004 

1 .98 

0.00 

5-23 

004 

3.63 

0.34 

6.55 

004 

4.27 

1.53 

8.  12 

005 

14.86 

4.41 

44.30 

005 

23.36 

6 . 06 

58.68 

005 

23.62 

7.66 

52.91 

005 

18.84 

4.72 

32.26 

005 

13.64 

0.80 

30.29 

005 

16.27 

1 .05 

44.60 

005 
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30  May  1975  (continued) 


Mean 

Maximum  No.  of 

Rate  Points 


T i me 


Ra i nfa 11  Rate 
(mm/hr) 


Mini  mum 
Rate 


1551 

1552 

1553 

1554 

1555 

1556 

1557 

1558 

1559 
1600 
160  1 
1602 

1603 

1604 

1605 

1606 

1607 

1608 

1609 

1610 
161  1 
1612 

1613 

1614 

1615 

1616 

1617 

1618 


22.09 

22.68 

16.94 

11.20 

10.40 

15.77 

11.74 

12.18 

9.82 

5.47 

5.04 
10.87 
12.23 

8.52 
6.19 
4.79 
4.57 
7-83 
6.14 
15.44 
13.28 
3-49 
8.  19 

8.05 
7.88 
7.  16 
7-17 
6.90 


6.75 

36.96 

005 

11.27 

38.10 

005 

4.83 

42.48 

005 

4.69 

23-02 

005 

1.11 

18.66 

005 

1.43 

37.29 

005 

1.88 

32.86 

005 

3.29 

35-46 

005 

2.38 

27.85 

005 

0.89 

8.28 

005 

0.91 

12.73 

005 

0.91 

42.35 

005 

1.38 

33.24 

005 

1.79 

20.93 

005 

1.51 

1 1.57 

005 

1.51 

8.03 

005 

0.96 

1 1 .87 

005 

0.87 

19-07 

005 

4.52 

8.06 

005 

3-52 

48.69 

005 

2.40 

49  .08 

005 

0.77 

8.98 

005 

2.91 

16.17 

005 

3-29 

11.57 

005 

3-07 

15.45 

005 

1.09 

15.88 

005 

1.14 

16.60 

005 

1.48 

20.03 

005 

NUMBER  OF  DROPS  PER  CUBIC  METER 


1 
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APPENDIX  D 

30  APRIL  1975  CASE  STUDY 

The  CHILL  radar  data  from  the  storm  of  30  April  1975  were  analyzed  for 
reflectivity,  liquid  water  content,  and  rainfall  rate.  Ground  data  were 
obtained  from  six  raingages  set  out  along  the  271°  radial  (see  Fig.  1 in  the 
report).  Liquid  water  content  was  determined  in  two  ways,  by  the  Ecc 1 es-Mue  I 1 e r 
attenuation  equation  (1971),  and  by  using  the  10  cm  reflectivity  in  a Z-M 
relationship  based  on  the  Jones  (1956)  central  Illinois  thunderstorm  data  and 
a relationship  given  by  Atlas  (1964).  Z and  M are  determined  in  the  following 
manner : 

M = 0.52R0-97 

1 4 1 

Z = 31500m1 

The  Jones  Z-M  equation  was  used  because  the  radar  data  discussed  here  were 
gathered  in  central  Illinois.  However,  a comparison  of  this  equation  with  that 
of  Marsha  1 1 -Pa I mer  and  India  thunderstorm  data  (Battan,  1973)  indicates  no-signi- 
ficant  difference  between  the  liquid  water  content  amounts  calculated  by  each 
method.  Table  1 contains  the  comparative  data  calculated  for  measurements  made 
at  gage  65  in  the  East-Central  Illinois  network  on  30  April  1975.  The  rainfall 
rates  are  5 minute  averages  centered  on  the  indicated  time. 

The  reflectivity  data  required  considerable  smoothinq  prior  to  use  in  any 
of  the  radar  equations  because  of  equipment  eccentricities,  and  the  real-time 
variability  of  the  precipitation  (that  is,  hydrometeor  scatterers).  Smoothing 
of  the  observed  reflectivity  data  for  the  calculation  of  liquid  water  content 
was  achieved  in  the  following  way. 

Three-bin  moving  averages  were  calculated  along  each  radial  between  269 
and  773(  azimuth  and  at  a given  elevation  angle.  The  averaged  reflectivity 
measurements  thus  generated  were  further  averaqed  by  accumulation  of  values  in 
each  bin  in  the  horizontal  plane  between  the  azimuth  constraints  and  in  the 
vertical  between  various  elevation  angle  intervals.  The  elevation  angle  intervals 
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used  were  0.3-0. 7,  1. 3-1.7,  2. 3-2.7,  and  3 - 3- 3 - 7 degrees.  These  sums  of  3_bin 
averaged  reflectivity  were  divided  by  the  number  of  measurements  yielding  the 

j 

final  values.  The  resulting  spatially  averaged  reflectivity  values  were 
attributed  to  a ray  path  centered  between  the  specified  azimuths  and  elevation 
angle  intervals. 

Radar  reflectivity  data  from  a 60  minute  period  between  1600  and  1700  CDT  on 
30  April,  1975  were  used  to  calculate  liquid  water  content  by  the  attenuation 
and  reflectivity  methods  and  vertical  cross-sections  were  plotted  and  analyzed. 

The  rain  system  moved  from  approximately  268°  nearly  parallel  to  the  271°  radar 
radial.  This  direction  of  movement  resulted  in  relative  motion  toward  the  radar 
with  little  cross-beam  speed.  The  motion  of  a raincell  associated  with  this  storm 
can  be  seen  in  Figures  1 to  11. 

The  liquid  water  content  was  calculated  by  the  Ecc les-Hue 1 ler  method  along 
o 

the  271  azimuth  by  using  the  range  difference  of  reflectivity  values  between 
successive  bins  to  calculate  attenuation.  The  final  liquid  water  content  estimate 
was  obtained  from  an  average  of  20  bin  values  centered  above  a gage  location. 

Similarily,  the  Z-M  equation  was  solved  for  liquid  water  content  in  each  bin 
followed  by  the  20  bin  averaging  centered  on  the  surface  gage  s i te . The  resu 1 ts 
are  discussed  more  fully  in  the  next  section. 

Liquid  Water  Content 

Figures  1 through  11  depict  liquid  water  fields  calculated  by  the  Z-M  relation- 
ship and  by  the  attenuation  as  calculated  from  the  dua I -wa ve 1 eng t h radar  at 
three  minute  intervals  over  one  hour's  time.  The  trianqles  in  the  figures  indicate 
the  location  of  the  raingages  along  the  271°  radial,  with  gage  71  closest  to  the 


L 


radar  and  gage  64  furthest  away.  The  dots  in  the  figure  indicate  the  points 
above  the  gages  where  the  final  estimates  of  liquid  water  content  were  computed. 
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Liquid  water  content  f i e I ds de term!  ned  by  using  the  attenuation  equation 
had  stronger  gradients  throughout  the  60  minute  period.  The  values  of  liquid 
water  ranged  from  less  than  1 g m~3  to  in  excess  of  9 g m The  liquid  water 
estimated  in  this  manner  seems  unrealistic  in  light  of  the  large  amount  of 
calculated  water  that  apparently  does  not  reach  the  ground.  For  example,  there 

_ O 

is  a closed  isoline  of  3-0  g m 3 above  gage  65  (between  65  and  70  km)  in  Figure  Id. 
This  area  of  high  liquid  water  remains  aloft  through  1 6 1 8 COT  (Figure  4b),  and 
there  is  no  evidence  from  the  raingage  data  of  this  amount  of  water  reaching  the 
ground  during  or  after  this  period.  At  1621  COT  abnormally  strong  liquid  water 
gradients  became  evident  betweeen  30  and  40  km  from  the  radar,  continuing  through 
1654  CDT  (Figure  10b). 

Continuity  between  liquid  water  calculated  aloft  and  ground  values  appears 
to  be  maintained,  however,  in  the  fields  calculated  by  the  Z-M  relationship.  The 
liquid  water  shown  in  Figure  1 c at  about  800  m above  the  second  gage  from  the  right 
appears  on  the  ground  3 minutes  later  (Figure  2a).  This  continuity  is  generally 
maintained  throughout  the  remainder  of  the  60  minute  period,  although  the  magnitude 
of  the  values  decrease.  There  is  also  a lag  ev i den t between  surface  and  above 
ground  values  (Figures  3a  and  3c).  This  may  be  caused  by  the  use  of  5~minute 
averaged  rainfall  rates  to  compute  the  liquid  water  at  the  gages  whereas  the  values 
above  ground,  derived  from  the  radar  reflectivity,  are  instantaneous. 

In  general,  the  liquid  water  calculated  from  the  rainfall  in  the  gages 
does  not  correlate  well  with  the  liquid  water  calculated  from  the  attenuation  of 
the  radar  beam.  The  radar  estimates  by  the  attenuation  technique  are  approximately 
five  to  six  times  higher  than  what  would  normally  be  expected. 

Liquid  water  content  calculated  using  the  Z-M  relationship  has  a much  better 


correlation  to  the  liquid  water  calculated  rrom  the  raingage  data.  There  were 
some  high  values  along  the  ground,  but  these  can  be  probably  be  accounted  for  by  the 
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the  five  minute  rainfall  rates  used.  The  gradients  of  liquid  water  by  the  Z-M 
relation  were  more  realistic  than  those  calculated  using  the  attenuation. 
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Marshal  1 


Jones  , Palmer  India-TSTM 


T i me 

Rainfall  Rate 

Z (dbz) 

M (9/ni  ) 

Z 

h 

Z 

M 

1545 

5.48 

36.98 

0.27 

34.83 

0.32 

33-82 

0.29 

1550 

23-70 

45.70 

1 . 12 

45.0 

1 . 16 

42.78 

0.97 

1555 

57.49 

50.97 

2.64 

51  . 16 

2.54 

48.2  1 

2.02 

1500 

58.36 

51.06 

2.69 

51.27 

2.57 

48.31 

2.04 

1505 

3.83 

34.85 

0.19 

32.34 

0.23 

31  .63 

0.2  1 

1610 

0.0 

0 

0 

0 

0 

0 

0 

1615 

0.0 

0 

0 

0 

0 

0 

0 

1620 

0.0 

0 

0 

0 

0 

0 

0 

1625 

0.16 

15.96 

0.09 

10.27 

0.01 

12.18 

0.02 

1630 

0.78 

25.38 

0.04 

21.28 

0.06 

21.88 

0.06 

1635 

0.74 

25-07 

.04 

20.92 

0.05 

21.56 

0.05 

1640 

0.54 

23.2 

0.03 

18.72 

0.04 

19.63 

0.04 

Table  I.  Z and  M data  calculated  by  three 
different  methods  from  rainfall  rate  data 
at  gage  65,  30  April,  1975- 
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Figure  1.  Liquid  water  content  fields  calculated  by  a Z-M  relationship 
(a  and  c)  and  by  radar  attenuation  (b  and  d)  for  1600:28  and  1603:26  CDT 
on  30  April  1975.  In  a and  b,  the  calculated  values  are  shown,  but  in  all 
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Figure  2.  Same  as  Figure  1 for  1606:20  and  1609:11  CDT 
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Figure  4.  Same  as  Figure  1 for  1618:02  and  1621:04  CDT 
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Figure  6.  Same  as  Figure  1 for  1630:03  and  1632:45  CDT 
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Figure  7.  Same  as  Figure  1 for  1635:51  and  1639:04  COT 
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Figure  8.  Same  as  Figure  1 for  1642:17  and  1645:30  CDT 
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Figure  9. 


Same  as  Figure  1 for  1648:24  and  1661:03  CDT 
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Same  as  Figure  1 for  1654:03  and  1657:06  CDT 
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